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Abstract 

Exclusive J /ip production, 7*p — > J /ipp^ offers a unique opportunity to determine the 
gluon density of the proton in the small x domain. We use the available HERA data to 
determine the gluon distribution in the region 10^^ 1z ^ ^ 10~^ and 2 ^ ^ 10 GeV^, 
where the uncertainty on the gluon extracted from the global parton analyses is large. 
The gluon density is found to be approximately flat at the lower scale; it is compared with 
those of recent global analyses. 



1 Introduction 

Global analyses do not reliably determine the gluon for x < a few 10~^ at low, yet perturbative, 
as shown in Fig. [1] This is due partly to the lack of precise structure function data for 
X < 10~^ and mainly due to the fact that the data included in global analyses actually probe 
the quark distribution, while the gluon density is constrained by the logQ^ dependence of the 
data, that is by the evolution. In the low x region the available interval decreases and the 
accuracy of the gluon determination becomes worse. The strong dependence of the global fits 
for the gluon on the order of the analysis is clearly demonstrated in Fig. [H Note that the recent 




Figure 1: Comparison of recent global fits of the gluon distribution at small x at leading 
(LO), next-to- leading (NLO) and next-to-next-to leading (NNLO) order, for the two scales 
Q'^ = 2.4 (left) and 10 GeV^ (right panel). LO gluons (dash dot) compared are CTEQ6L [1] and 
MRST2004F4LO [2]. The two (long) dashed lines indicate the error estimate of the CTEQ6.5 [3] 
gluon and the shaded band is the error band for the MRST2001 [1] global gluon. Central values 
for the NLO global fits are from CTEQ6.5M (short dashed) and MRST2004NLO [5] (dotted). 
The solid line represents MRST2006NNLO [6]. 
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gluon from the MRST NNLO analysis [6] receives sizeable corrections both in size and shape 
compared to the NLO fit, signalling a large uncertainty of the gluon in this regime. In this 
context it is also interesting to note that the gluon as obtained from global fits can significantly 
change, both in normalisation and shape, if small x resummations are incorporated into the 
analysis [7]. 

Data for the exclusive 7*p —^J/ipp process offer an attractive opportunity to determine the 
low X gluon density in this domain, since here the gluon couples directly to the charm quark 
and the cross section is proportional to the gluon density squared [S]. Therefore the data are 
much more sensitive to the behaviour of the gluon. The mass of the cc vector meson introduces 
a relatively large scale, amenable to the perturbative QCD (pQCD) description not only of 
large diffractive electroproduction, but also photoproduction of J/ip. The available J/ip 
data probe the gluon at a scale /i^ in the range 2 — 10 GeV^ for x in the range 10~^ ^ 2; < 10^^; 
that is just the domain where other data do not constrain the gluon reliably, see Fig. [H It 
would be good to have comparable data on exclusive T production to determine the gluon at 
larger scales, but here the available data are sparse, see Fig. [5] below. 



2 Exclusive J/i/j production at LO 

To lowest order the 7*p ^ J/tpp amplitude can be factored into the product of the 7 cc 
transition, the scattering of the cc system on the proton via (colourless) two-gluon exchange, 
and finally the formation of the J/ ip from the outgoing cc pair. The crucial observation is that 
at high 7p centre-of-mass energy, W, the scattering on the proton occurs over a much shorter 
timescale than the 7 ^ cc fluctuation or the J/ ip formation times, see Fig. [2l Moreover, at 
leading logarithmic accuracy, this two-gluon exchange amplitude can be shown to be directly 
proportional to the gluon density xgi^x^Q"^) with 



= (Q2 + Mj/^)/4, 



X = {Q' + M'j/^)/{W' + Mj/^] 



(1) 



is the virtuality of the photon and Mj/^ is the rest mass of the J/ip. To be explicit, the 
lowest-order formula is ^ 



J/ip p) 



t=o 



48a 



-xg{x,Q^ 



1 + 



(2) 



where Fee is the electronic width of the J/ip. 

In the leading logarithmic approximation, the integral over the transverse momentum kx of 
the t-channel gluons, see Fig. [21 gives rise to the integrated gluon density g{x, Q"^). As usual in 
collinear factorisation, the kx dependence of the integral is completely absorbed in the input 
gluon distribution (of the global analyses), taken at the factorisation scale Q^. The integral 
over the charm quark loop is expressed in terms of the electronic width. Fee, of J/^/', and 
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Figure 2: Schematic picture of high energy elastic J /ip production, 7*p — > J /ipp. The factorised 
form follows since, in the proton rest frame, the formation time Tf ~ 'iE^/ (Q^ + Mj^^) is much 
greater than the cc-proton interaction time Tint- In the case of the simple two gluon exchange 
shown here. Tint — -R, where R is the radius of the proton. 

the Q^/Mjy^ term in the final brackets reflects the contribution of the longitudinally polarised 
incoming 7*. Equation ([2]) gives the differential cross section at zero momentum transfer, t = 0. 
To describe data integrated over t, the integration is carried out assuming a ~ exp(— 6t) with 
h the experimentally measured slope parameter. Throughout this work the value 

b = 4.5 GeV-2 (3) 

is used, which is in ag reement with P [101 EjQ Thus it becomes possible to extract the gluon 
density g{x,Q'^) directly from the measured diffractive J /if) cross section. 

3 Corrections to LO exclusive J/ip production 

Expression ([2]) is a simple, first approximation, justified in the LO collinear approximation using 
the non-relativistic J/ip wave function. The relativistic corrections were intensively discussed 
in [T2I |13]. The problem is that, simultaneously with the relativistic description of the c 
quarks, one needs to consider higher order Fock component ccg states of J/ip. Hoodbhoy [H] 
has studied these two effects to order f^/c^. He has shown that relativistic corrections to ([2]), 
written in terms of the experimentally measured Fge, are small, ~ (9(4%), see [Hj. So we do 
not account further for relativistic corrections below. 

NLO corrections arise, first, from an explicit integration over the gluon kx-, which goes 
beyond the leading log contribution, arising from dk^/k'^, and, second, from more complicated 

^We neglect a slight energy dependence of b which is only observed for photoproduction, but which is of the 
order of differences between measurements of the two experiments HI and ZEUS. The possible uncertainty is 
smaller than, or at most comparable to, other approximations used. 
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diagrams with one additional loop. To perform the explicit kx integration we have to use the 
unintegrated gluon distribution, /(x, which is related to the integrated gluon by 



xg{x,fi')= I ^/(x,4) + c(Q2). (4) 







T 



Of course, the infrared contribution cannot be treated perturbatively, and so we have introduced 
a lower limit Qq for the integral in the J/ip production amplitude. 



g4 xg{x,Q 



2^ T HH^^-Mj/^)/4 dA;2«,(fc2) d\xg[x,k\)\ , aM 



Q\Q^ + kl) dkl + g4 ^9{x.Q,). (5) 



Expression ([5]) replaces the factor as{Q'^)xg{x,Q'^)/Q^ in the LO result ([2]). To be precise, the 
unintegrated distribution / embodies the Sudakov factor T(/c|,,/i^) such that 



fix, kl) = d [xg{x, kl)T{kl, f?)] Id In k% . (6) 

Thus c((5o) in Eq- Q is given by xg{x,Q'^T{Q\,ix^) and correspondingly for ([5]), xg{x,Q'^ — >• 
Qg)T(QQ, yu^). In our numerics we have chosen [x^ = Q"^. However, in the amplitude 
([5]), the dominant contribution comes from the region oi kx ^ Q where T(fc|.,/i^) is close to 
unity. The inclusion of the T factor may be considered as an 0{as) correction to the gluon 
density and suppresses the gluon in our analysis by 1.7% for photoproduction at x = 10~^. The 
contribution coming from kx < Qo is written in terms of the integrated gluon g{x, Ql), that is 
the infrared part is absorbed into the input distribution at the 'transition' scale Qo- 

Of course, at low the gluon extracted from a global analysis may be affected by the 
presence of absorptive corrections which are usually neglected. Here, the absorptive corrections 
are expected to be smaller. The transverse size, r, of the qq dipole produced by the 'heavy' 
photon in DIS has a logarithmic distribution J dr^/r^ starting from up to some hadronic 

scale. In the case of J/ip production the size of the cc dipole is limited by the size of the J/ip 
meson. Even in photoproduction it is of the order of Since the probability of rescattering 

is proportional to r^, we anticipate a much smaller absorptive effect. 

A more detailed analysis of the NLO corrections was done in fi6[ [T7] . Part of these cor- 
rections generates the running of a^, while part is similar to gluon Reggeization in the BFKL 
approach. Indeed, for J/ip electroproduction it was shown [TTj using the conventional coUinear 
factorisation scheme, that there is a NLO correction of the form 

In the kx factorisation approach such a correction may be included by replacing the t-channel 
gluon by the Reggeized gluon. However this correction vanishes with a natural choice of the 
factorisation scale, = Q^, which was adopted in our prescription [HI [19], [15]. One therefore 
has reason to believe that the kx factorisation approach accounts for a major part of the NLO 
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effect, and that the resulting 'NLO' gluon may be compared to that in a set of NLO global 
partonsl^ Therefore we shall refer to the resulting distributions as NLO gluons. 

One also needs to account for the fact that the two gluons exchanged carry different fractions 
X, x' of the light-cone proton momentum, see Fig. [21 That is one has to use the generalised 
(skewed) gluon distributionjfl In our case a;' <^ a; <^ 1, and the skewing effect can be well 
estimated from [201 

22^+3 r(A + |) 

where A((5^) = d [ln(x(7)] /91nx. That is in the small x region of interest we take the gluon to 
have the form xg ~ x^^ . 

Recall that the integral ([5]) was written for the discontinuity (i.e. for the imaginary part) of 
the amplitude shown in Fig. [2J The real part may be determined using a dispersion relation. 
In the low x region, for our positive-signature amplitude A oc x~'^ + (— x)~^, the dispersion 
relation can be written in the form 

ReA TT, TT SlnA TT^ln (x5f(x, Q^)) 

~ — \ ^ ^ L (Q) 



ImA 2 2 91n(l/x) 2 d\n{l/x 
Both corrections lead to an enhancement of the cross section 



4 Determination of the gluon from J/i/j data 

In the following, we present fits to the data for exclusive J/ ip production from HERA using the 
perturbative description discussed above. 

In the low x region it is expected that the x dependence of the gluon density xg{x, Q"^) is well 
approximated by the form x~^. However, the evolution in Q"^ modifies this behaviour, enlarging 
the power A as increases. In particular, in the double leading log (DLL) approximation, we 
have the asymptotic form 



xg ~ exp ( \j ln(l/x) In Q2 j . (lo) 



^The global partons are defined in the MS regularization scheme. Our partons should also be considered to 
be in the MS scheme, since we use the MS definition of a^, and moreover the factorisation scale which provides 
the cancellation of the as\n.l/x correction is also specified in the MS scheme. 

^In the formal analysis of the NLO contributions, there are effects arising from integrated quarks which 
generate gluons which then couple to the charm quark. In terms of our unintegrated gluon (/) description this 
should be considered as a NLO correction to the evolution of /. Since we do not consider the evolution, but 
just parametrise the scale dependence of the gluons (sec below), this correction is outside our analysis. 
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Thus we need a dependent parametrisation. However, in the hmited region of covered 
by the exclusive J/ip data, it is sufficient to use a simple parametric tornfl 



-X 



with X = a + b ln(/iV0.45 GeV^) . 



(11) 



The free parameters N, a and b are determined by a non-hnear fit to the exclusive J/ip 
data from HI [9] and ZEUS [lT],[22]if| This three-parameter form provides enough flexibility to 
accurately describe the x and behaviour of J/ip production in the limited domain covered 
by the J/ijj data, namely 10"^ < a; < lO'^ and 2 < < 3 QeV^, so we will use exactly 
Eq. ([n]) for the LO fit. 

However, for the NLO approach, where we have the fc^ integral ([5]) which runs up to 
the kinematical limit = {W"^ — Mj^,^)/A, we face the problem of a badly convergent 
integral. Indeed the low x gluons xg ~ obtain a large anomalous dimension 7. Since 
^-bink^ _ ^^2^^f)in(i/a;) ^ _ JjIy^I^i/x). On the other hand, the expression was calculated 

with running a^. Recall that the scale dependence of the power A was generated by the evolution 
of the form A = J as{q'^)^- Accounting for the running a^, it is natural to replace the second 
term in the parametric form (fTTl) by Inlnlfi"^ / Aqq^^) . Therefore, for the NLO fit, we use 



The gluon densities obtained from the LO fit, using ([2]) with ( |TT|) . together with the skewing 
and real part corrections, to the exclusive J/ ip data [H |TTl are shown in Fig. [3]@ Inclusion of 
these corrections gives a 22% suppression of the gluon for photoproduction at x = 10"'^, with 
the skewing correction giving the dominant suppression of 18%. We use a 1-loop running 
with as(M|) =0.118. In the analysis, error bands on the gluon and cross section are generated 
using the full covariance matrix for the fitted parameters, where as input we have added the 
statistical and systematic experimental errors of the data in quadrature. Compared to the 
gluons from the global fits, the gluon from our LO fit is similar in shape, but slightly smaller 
in normalisation and shows less rise towards smaller x with growing scales. 

We now present the results of our NLO fit, which we obtain by modifying Eq. ([2]) by help 
of ([5]) with 2-loop running a^, and replacing the parametrisation ffTTl) by f|T2l) . Of course we 
also include skewing and real part corrections as before, and the T factor as discussed above. 
Figure H] shows our fit using Ql = 2 GeV^ and A^cd = 0-09 GeV^ The dotted lines represent 
the central values of the gluons obtained using Qg = 1 GeV^. The NLO gluon fit shows a 
better matching to the global gluons at x = 10^^ than the LO gluon obtained. The analysis 

''Such a form has already successfully been used in [21] for the analysis of inclusive diffractive DIS data. 

^We also performed fits for the data from the HI and ZEUS collaborations separately. These fits typically 
have a smaller X^iu, signalling a slight incompatibility between the data. However, they lead to similar results 
for the gluon. As the combined fit is very satisfactory (x^j^^/d.o.f. < 1), we will not discuss fits of individual 
data sets in the following. 

^Note that in Figs. [21 [H only a subset of 51 data points used in the fits is displayed, and there are data points 
at up to (Q^) = 22.4 GeV^, corresponding to = 8 GeV^. 



xg{x,fi'^) = Nx 



-A 



with \ = a + b Inln^fi^ / Aqqj)) . 



(12) 
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Figure 3: Leading order fit of elastic J /iIj data as described in tlie text. Left panel: cross section 
compared to some of the HI p] and ZEUS [HI [22] data, with values for as indicated; right 
panel: gluon compared to global fits for scales as indicated, where the solid (dashed) lines are 
the MRST2004F4LO [2] (CTEQ6L [1]) results. The width of the bands displays the uncertainty 
of the cross section and fitted gluon respectively, whereas the darker shaded areas indicate the 
region of the available data. 







a 


h 


Xmin/d-o.f. 


LO 


0.99 ±0.09 


0.051 ±0.012 


0.088 ±0.005 


0.9 


NLO 


1.55 ±0.18 


-0.50 ±0.06 


0.46 ±0.03 


0.8 



Table 1: Values of the three parameters of the LO and NLO gluon fits and corresponding 
Xmin/d-o.f. 
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W (GeV) X 

Figure 4: As Fig. [3], but for the next-to-leading order fit and comparing to NLO global fits from 
CTEQ6.5M [3] and MRST2004 [5]. 



Q2(GeV2) 




Amrst 


-^CTEQ 


2.4 


0.04 


-0.17 {-1.07, -0.16, 0.00} 


0.01 {0.04, 0.00, 0.05} 


4.1 


0.11 


0.06 {-0.03, 0.07, 0.16} 


0.13 {0.14, 0.13, 0.19} 


6.4 


0.16 


0.15 {0.09, 0.15, 0.24} 


0.19 {0.18, 0.19, 0.27} 


8.0 


0.19 


0.18 {0.13, 0.18, 0.27} 


0.21 {0.20, 0.22, 0.30} 



Table 2: The values of the power of the gluon. A, at four values, for our NLO fit to 
elastic J/ production data compared to two global fits [U [5] . The numbers for MRST and 
CTEQ are obtained through a fit in the range x = 10~^ . . . 10~^, assuming xg ~ with an 
X independent A, whereas the values in curly brackets are the logarithmic derivatives of the 
gluons at X = {10~^, 10~^, 10~^}, respectively. 
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of the exclusive J /if) data indicates that, at the larger scales, the small x behaviour of the 
gluon distribution is slightly flatter than that of the global analyses, both in their x behaviour 
and in their scale dependence of A(Q^). At low scales, the gluon obtained from the J /if) data 
still rises with decreasing x, especially in contrast to the MRST fit. For completeness we also 
present the values of the parameters for the LO and NLO fits (Tab. [1]). As can be seen from 
Figs. [31 [Hand the Xmin/^-o.f. values quoted in Tab. [T], our simple ansatz for the form of the 
gluon, xg ~ using an x independent power A, works very well. To quantify the x and 
behaviour, we tabulate the values of the power of the gluon. A, from (fT2|) at four values, 
compared to values estimated from MRST2004NLO [5] and CTEQ6.5M gluons [3] (Tab. W)- 
As is evident from Fig. H] and Tab. [21 our NLO gluon fit seems to be incompatible with the 
strength of evolution of the MRST global fit (e.g. at x = 10"'^ our gluon increases by a factor 
2.7 from = 2.4 to 8 GeV^ compared to 3.7 for the MRST2004NLO gluon); however, there 
is fair agreement in the evolution between our gluon and the CTEQ fit (CTEQ evolves by a 
factor 2.8 in the same regime, although in absolute normalisation our NLO J /if) prediction for 
the gluon agrees, on average, much better with the MRST NLO prediction). Of course, this 
should be seen in light of the large uncertainties of both the MRST and CTEQ gluons at small 
X and scales, see Figs. [T]and[ll This uncertainty persists at the largest scales probed in our fit. 
As these scales are rather low, the accuracy of the DGLAP approach may already be seriously 
affected by small x effects and absorptive and power suppressed corrections. 

Of course, in the kx factorisation approach there is some uncertainty arising from the 
infrared cut-off Qo? below which we cannot consider the kx integration literally. We have to 
express this low kx contribution in terms of the gluon integrated over k"^ up to Qq. However, 
with our prescription for unintegrated partons, these two contributions match smoothly to each 
other. The ambiguity due to the choice of Qo is quite small, as illustrated by the closeness of 
the dotted lines to our NLO gluons in Fig. [H 

Note that the difference between the LO and NLO gluons is large at the smallest x values, 
both in the global parton analyses and for the gluons obtained from elastic J/ip production. 
For the global analyses this is due to the absence of the photon-gluon coefficient function at 
LO. At LO, the photon couples only to the quark parton, which is produced from a gluon 
at some scale <^ Q^, due to the strong ordering in kx- At low x the gluon grows with 
g^, and, to provide the measured values of the proton structure function F2, we need a much 
larger gluon distribution in the LO formalism. At NLO the photon-gluon coefficient function 
is present, which provides a direct photon-gluon parton coupling at the scale Q^, and, more 
importantly, a. 1/x divergence appears in the quark-gluon splitting function, which accelerates 
the quark evolution and in turn requires less gluon. In elastic J/ip production we have an 
analogous situation. By carrying out the kx integration, we include the interaction with the 
gluon at large scales of the order of + Mj^^. Moreover, part of this integral has k^ ^ Q^, 
and may be regarded as one step of backward evolution. In summary, in DGLAP analyses 
based on collinear factorisation, the large change in the gluon distribution in going from LO to 
NLO is due to the strong ordering in kx and the absence of an additional loop integral in the 
LO coefficient function and parton evolution. Inclusion of higher-order terms beyond NLO is 
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Figure 5: Prediction of elastic T photoproduction, using our LO and NLO gluon, compared to 
data. The dotted lines indicate the error band of the LO prediction, whereas the shaded band 
is our NLO prediction. The data points are the published ZEUS [23] and HI [21] results, and 
the preliminary results from ZEUS [25] . 



expected to give a much smaller effect. They will mainly affect the normalisation and not the 
X dependence of the gluon. 

The higher-order corrections in our approach are 0{as), but may be enhanced by large 
logarithms. The corrections which are enhanced by ln(l/x), and which may lead to some x 
dependence, are absorbed in the form of the gluon distribution by choosing the appropriate 
factorisation scale //^ = Q^, see Eq. ([7]). Then, in our NLO approach, the part of the 
integral which may be logarithmically large, is accounted for explicitly in the factorisation 
formalism which is used to obtain the fcy integral. After this, those higher-order corrections, 
which are not included in the /ct integral, are concentrated in the domain fc^ ~ A* = <5- Thus 
the scale dependence of these higher-order corrections is driven mainly by the scale dependence 
of the running of which we take into account. The fact that the gluons obtained from our 
analysis turn out to be close to the gluon distributions coming from the global analyses for 
X ^ 10~^, where they are fairly stable, indicates that the omitted higher-order corrections are 
indeed small. 

We have extended our framework to predict elastic T photoproduction, using our LO and 
NLO gluon with our cross section formulae including corrections and T factor, see Fig. O 
Although the data is sparse, the cross section predictions are reasonable. 

We conclude that this new information coming from our NLO analysis of elastic J/ip pro- 
duction data, in which the cc couple directly to the low x gluon parton and where the cross 
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section is proportional to the square of the gluon, is especially valuable to constrain the small 
X behaviour of the gluon distribution. The accuracy of the elastic J/ip data is now sufficient, 
as indicated by the error bands (arising from the experimental uncertainties) on the extracted 
gluon distribution shown in Fig. HI to improve our knowledge of the gluon distribution at small 
X, X ^ 10~^, considerably. This is in comparison to the small x behaviour of the 'global' gluon 
distributions, which is not well constrained by the inclusive DIS structure function data. 
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